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Abstract 

This paper proposes a novel nature-inspired optimization algorithm called the Fox optimizer (FOX) which 
mimics the foraging behavior of foxes in nature when hunting prey. The algorithm is based on techniques 
for measuring the distance between the fox and its prey to execute an efficient jump. After presenting the 
mathematical models and the algorithm of FOX, five classical benchmark functions and CEC2019 
benchmark test functions are used to evaluate its performance. The FOX algorithm is also compared against 
the Dragonfly optimization Algorithm (DA), Particle Swarm Optimization (PSO), Fitness Dependent 
Optimizer (FDO), Grey Wolf Optimization (GWO), Whale Optimization Algorithm (WOA), Chimp 
Optimization Algorithm (ChOA), Butterfly Optimization Algorithm (BOA) and Genetic Algorithm (GA). 
The results indicate that FOX outperforms the above-mentioned algorithms. Subsequently, the Wilcoxon 
rank-sum test is used to ensure that FOX is better than the comparative algorithms in s statistically 
significant manner. Additionally, parameter sensitivity analysis is conducted to show different exploratory 
and exploitative behaviors in FOX. The paper also employs FOX to solve engineering problems, such as 
pressure vessel design, and it is also used to solve electrical power generation: economic load dispatch 
problems. The FOX has achieved better results in terms of optimizing the problems against GWO, PSO, 
WOA, and FDO. 
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1. Introduction 

Currently, meta-heuristic optimization algorithms are widely used to solve real-world application problems [1]. These 
algorithms are predominantly inspired by biological or physical behavior. Algorithms that are inspired by physical 
behavior include space gravitational optimization (SGO) [2] and the gravitational search algorithm (GSA) [3]. These 
algorithms are derived from the theory of relativity, the law of gravity, and interactions between objects. 


Consequently, different metaheuristic algorithms have been proposed by researchers since 1948, when Alan Turing 
broke the code of the Enigma cipher machine. After the idea of Turing’s heuristic method, the genetic algorithm 
(GA)[4] was proposed, which mimics natural evolution. After the GA was proposed, many algorithms have been 
developed, such as simulated annealing (SA) [5], tabu search [6], ant colony optimization (ACO) [7], particle swarm 
optimization (PSO) [8], and differential evolution (DE) [9]. Some of the recent algorithms in this family are dragonfly 
optimization algorithm (DA) [10], whale optimization algorithm (WOA) [1], grey wolf optimization (GWO) [11], 
butterfly optimization algorithm (BOA) [12], cat swarm optimization (CSO) [13], the fitness-dependent optimizer 
(FDO) [14], [15], the donkey and smuggler optimization algorithm (DSO) [16] and the learner performance-based 
behavior algorithm (LPB) [17], K-means GWO [18], chimp optimization algorithm (ChOA) [19]. Each algorithm has 
different inspirations and operators. For example, GA uses three basic operators: selection, crossover, and mutation 
[20]. The GA and its operators are used by researchers to solve optimization problems. Genetic algorithm is quite 
efficient but the problem presentation is always a challenge for researchers due to the use of discrete chromosomes. 
GA also has premature convergence and it has many parameters to set, such as population size, crossover rate, 
mutation rate, and the selection method [21][22]. 


PSO [23] was proposed by Kennedy in 1995 that mimics the behavior of swarm particles. Each particle represents a 
solution in the search space [24]. PSO uses the global best position of the particle and the local best position of particles 
to calculate new velocity. The new velocity is based on the memory of the previous particles and the neighbor’s 
experience. Then, the new position of the particle is found based on the previous position and the velocity [25]. Despite 
using PSO for solving problems in different areas including neural network training, estimation, and face recognition, 
it has disadvantages of falling into local optima and having low convergence when dealing with high dimensional 
problems as discussed in [26], [27], [28]. 


DA is a meta-heuristic algorithm inspired by the static and dynamic swarming behavior of dragonflies in nature [29]. 
DA is using these behaviors either for migration or hunting. DA has three basic operators to explore and exploit food 
sources. These operators are separation, alignment, and cohesion [10]. DA had been used in different areas to solve 
real-world applications: image processing, machine learning, and networking. Results proved that DA outperformed 
against PSO for solving multidimensional functions. DA is also superior to PSO while it is evaluated on CEC2019 
benchmark functions. According to [10], DA achieved competitive results against common algorithms. DA also has 
drawbacks due to trapping into local optima and premature convergence. Using levy flight also makes more 
interruption in the search space due to the implementation of big steps. It also uses more time in the exploitation phase. 


After DA, another algorithm was proposed which is called WOA. Humpback whales are the source of inspiration for 
this algorithm proposed in 2016 [1]. WOA consists of different methods to catch prey. Firstly, it uses an encircling 
mechanism to create a search space around the prey. Then, it uses the bubble net attacking method to shrink the circle 
and update whales’ position toward the prey. The exploratory method is activated due to discovering better prey based 
on a random walk and best positions [30]. WOA had been used to solve many applications in different fields, such as 
electronic engineering, economic scheduling, civil engineering, fuel and energy, and medical engineering. Such 
diverse application areas demonstrated that WOA achieved high performance in terms of exploration and exploitation. 
WOA also has a high chance of escaping local optima by balancing both the exploration and exploitation phases. 
WOA has limitations regarding convergence speed and achieving optimum results. One of the problems is the 
randomization methods in WOA, these methods increase the computational time when it deals with a problem that 
has many dimensions. Having poor convergence is another disadvantage of WOA because both exploration and 
exploitation phases depend on a single variable. Besides, WOA tends to be trapped into local optima in large-scale 
problems because it uses an encircling mechanism[30]. 
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BOA is inspired by the butterflies foraging behavior. This behavior is responsible for moving the butterflies towards 
the food position. BOA was proposed in 2018 [12]. Butterflies discover the direction of the food source via their 
senses and analyze the smell of the food in the air. BOA consists of three phases which are the initialization phase, 
iteration phase, and final phase. The searching process is done in the iteration phase while the best position is found 
the final face is terminated [12]. BOA has been used to solve different applications: spring design, gear train design, 
and welded beam design problems. It is also used to solve the green vehicle routing problem. Results proved that BOA 
obtained high performance against PSO, Artificial Bee Colony (ABC), and Differential Evolution (DE) [12]. Despite 
having better performance, BOA has difficulties due to low convergence [31]. 


FDO is a recent meta-heuristic algorithm that was inspired by the reproduction process of bees while searching for 
better hives [14]. The main aim of FDO is how scout bees try to find better food sources among various hives. FDO 
works based on two parts which are the pace and positions of the scout bees. The pace is the velocity of the scout bees 
toward the food source. Therefore, the updating positions of the scout bees are based on the new velocity and the 
previous position of the scout [32]. FDO is a state-of-the-art algorithm used to solve different applications: antenna 
design, economic load dispatch problems, and power systems [33]. It achieved superior results compared to other 
common algorithms: GA, PSO, DA, ChOA, WOA. FDO has an O(n) time and space complexities. Despite achieving 
these results, there are still chances for improving FDO in terms of improving solutions [32], [15]. In addition, slow 
convergence and not having a proper balance between exploration and exploitation are disadvantages of FDO [34]. 
Therefore, it can be said that active nature-inspired algorithms entail being proposed by researchers to solve the issues 
of the state-of-the-art algorithm. 


Chimp Optimization Algorithm (ChOA) is a recent metaheuristic algorithm that mimics the behavior of chimps while 
they are hunting in a group proposed in 2020 [19]. It works based on individual intelligence and the sexual behavior 
of chimps. Therefore, four types of chimps are used in ChOA which are attacker, barrier, chaser, and driver. These 
types of chimps are responsible for exhibiting diverse intelligence. ChOA also has four main steps to hunt preys which 
are driving, chasing, blocking, and attacking [35]. ChOA has been used to solve many applications. According to 
[19] ChOA was used to solve ten design problems. ChOA achieved the best results against GA, PSO, and GWO for 
solving benchmark functions. However, it has limitations due to the transition between the exploration and 
exploitation phases [36]. It also has issues solving binary problems [37] and the solution requires more enhancement 
[35]. 


Lion Optimization Algorithm (LOA) was inspired by lion’s behavior and social interaction. LOA use social 
mechanism and encircling mechanism to hunt prey [38] while TIGER algorithm uses eight different technics to catch 
prey such as birth and breed, migration of children, migration, death, suicide, prowl, mating, and signing [39]. 


These metaheuristic algorithms mentioned above are used to optimize real applications in various fields, such as 
structural optimization [18], [40], scheduling and routing [41], software testing [42], image processing [43], and data 
mining [44]. However, metaheuristic algorithms may have issues related to their poor search capability and falling 
into local optima [45], [46]. 


For example, the WOA falls into local optima when confronting high-dimensionality problems because of having only 
an encircling mechanism in the exploration phase [47]. Therefore, this problem has been solved by using a new 
method, which is used to update the position and enhance the capability of the exploration phase [48]. In [47], chaotic 
local search is combined with the WOA to enhance the capacity of the WOA to overcome local optima. 


GWO is another example that falls into local optima when solving high-dimensional problems and complex 
multimodal problems [49], [50]. However, this problem is solved by using weight variables for different types of 
wolves in GWO. By using the weight variables and changing the control parameter a, GWO can achieve the maximum 
admissible error (MAE), which is required in engineering optimization problems [49]. 


The bat algorithm (BA) also has difficulties with some multidimensional functions because it can be trapped in local 
optima. Thus, this problem is solved by modifying the exploitation phase by changing the dimensional size to a 
dynamic dimension and adding an inertia weight [50]. 
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Additionally, having an improper balance between the exploration and exploitation phases is problematic in WOA 
[51], [52]. As a result, an improved WOA is proposed in [53] and hybridized with DE’s mutation operator. The 
modified WOA is called improved WOA (IWOA). The mutation operator is integrated inside the exploration phase 
of the WOA, and the exploitation phase of the IWOA is the same as that of the WOA. Besides, in [51], researchers 
solved this issue by hybridizing GWO with the WOA to improve the performance of the WOA and achieve a better 
balance between the two phases. Like WOA, GWO has the same problem. Consequently, researchers proposed 
improved GWO combined with DE to achieve a proper balance between the phases [54]. Therefore, the position of 
the wolves is updated based on the principle of the survival of the fittest (SOF). After adding the operation of DE, the 
position is updated in GWO. Similar to GWO, the BA also has this problem. Thus, an inertia weight is added to update 
the velocity of the bat to improve the efficiency of the BA in terms of balancing both phases [55], [51]. The inertia 
weight dynamically decreases throughout the iterations. It has a high value, in the beginning, to allow the BA to 
perform a global search. Decreasing the value allows the BA to have better local searchability. 


It should be noted that there are other fox-inspired algorithms in the literature. However, none of them used the solo 
hunting strategy modeled in this work. To the best of our knowledge, there are two papers published in the literature 
that attempted to propose an optimization algorithm inspired by fox's behaviours. One of them was published in 2016 
under the name fox hunting algorithm (FHA) [56] and the other was published in 2021, which is called the red fox 
optimization algorithm (RFO) [57]. FHA mimics the hunting behavior of the hunters who tries to hunt the red fox by 
using horses and hounds. Therefore, the FHA is related to the way the hunter wants to catch the fox. The RFO also 
uses the behavior of searching and hunting mechanism of red fox. However, RFO uses the reproduction mechanism 
and leaves behind the herd by red fox. RFO also uses a sorting mechanism inside the exploration phase while FOX's 
exploration is based on the best position and random search. Despite proposing both FHA and RFO, FOX is different 
from them because it uses the behavior of red fox for searching and hunting in the snow environment. It also mimics 
the jumping technic in order to dive into the snow to hunt the prey. In the exploitation phase, distance of red fox is 
measured by the prey and then it estimates the jump technic in order to dive and hunt the prey. FOX also requires an 
artificial fox to jump in northeast direction and its opposite. Therefore, the new position is found based on the distance 
of the prey, jump value, and the direction range value. 


The main aim of this paper is to overcome the stagnated local optima problem and not have the proper balance between 
exploration and exploitation phases. Therefore, the proposed FOX algorithm effectively balances the exploration and 
exploitation phases better than its alternatives. FOX is evaluated using CEC2019 benchmark functions. It can also be 
used to solve various optimization problems. Additionally, decreasing a specific parameter provides support for 
avoiding local optima. Using the exploration phase based on the best position and average time variable has a great 
impact on not trap into local optima. 


Therefore, the main contributions of this paper are: 


1. Proposing the FOX optimization algorithm based on fox hunting behavior. 

2. Solving the pressure vessel design problem and economic load dispatch problem. 

3. Using a unique technique in the exploration phase that can be adapted to propose new algorithms or modify 
algorithms. 


The rest of this paper is organized as follows: Section 2 describes fox life, and fox behavior, and then, the FOX is 
described in detail in Section 3. Section 3 describes the results and gives a discussion of the FOX, and then the pressure 
vessel design problem and economic load dispatch problem are explained and solved by the FOX. Finally, the 
conclusion and future work are presented. 


2. Inspirations, mathematical models, and the algorithm of FOX 

In this paper, a novel red fox in the snow optimization algorithm is proposed which is inspired by the hunting 
behavior of red fox. The algorithm uses both the searching and hunting techniques of the red fox. The following 
subsections present red fox life, red fox behavior in the snow, and the FOX algorithm. 
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2.1 Fox life 


A fox can live in a species-poor environment that has low production [58]. Foxes have different colors, but the most 
common colors are white and red. Two common types of foxes are the red fox and arctic fox. The red fox is the most 
widespread animal, and it has colonized urban areas in the USA, Europe, Canada, Japan, and Australia [39, 40]. The 
arctic fox has a strong ability to hunt its prey either from above or below. Fox diets include geese and lemmings in 
North America. However, a large part of the fox diet is lemmings in Fennoscandia, Siberia. The fox also lives on 
islands such as the Pribilof Island and Svalbard. On these islands, the fox focuses on resources in the sea, such as fish 
and birds [59]. According to [59], foxes that live on these islands have acquired a mouth that opens wider to catch 
more prey. The main difference between these two foxes is that Arctic foxes live at lower temperatures compared to 
red foxes. Arctic foxes have difficulties living in their habitat because of human interference [60]. 


The red fox is an example of a fox that eats mammals and plants. It lives between approximately 2 to 4 years, and the 
weight of this animal is between 6.5 and 24 pounds [61]. It lives in diverse locations, such as grassland, forest, desert, 
and mountains [62]. 


By being an omnivore, the fox can adapt itself to its habitat in terms of food resources. Consequently, the fox eats 
vegetables, fruit, frogs, and fish [63]. However, the arctic fox eats lemmings and geese [64]. Foxes use intelligent 
methods to hunt their prey [65]. 


2.2 Red fox behavior in snow 
The FOX algorithm mimics the hunting behavior of a red fox when it dives into the snow to hunt its prey. The main 
techniques are based on the red fox trying to hunt the best prey [64], [66]. The steps are described below: 


1. When snow covers the ground, the snow blocks the sight of the prey. Red fox tries to search for prey 
randomly. 

2. Consequently, the red fox finds the prey by hearing ultrasound from the prey. Then, it takes time to get to 
close the prey. 

3. By listening to the sound of the prey and the time difference, the red fox can determine the distance between 
itself and the prey. 

4. After determining the distance, the red fox estimates the jump that is needed to catch the prey. 

5. Walking is performed randomly according to the minimum time and best position. 


In the beginning, the red fox moves randomly in the search space to discover prey. It can find prey via hearing 
ultrasound of the prey. This random walk is what we inspired and use for providing exploratory behavior in FOX. 
While the red fox searches, it may hear the sound of the prey. After hearing the sound, the red fox is in the exploitation 
phase. The red fox can hear ultrasound, so the sound of the prey takes time to reach the red fox. The distance sound 
travels can be calculated by multiplying it by the speed of sound travel in the air, which is 343 meters per second. 
However, because of having static value, a different method is used to investigate the speed of sound travel in the air. 
The red fox tries to move forward the prey decides to jump against the prey. Therefore, the red fox tries to jump 
depending on the time the sound takes to catch prey. A study showed that a fox prefers to jump northeast based on the 
magnetic alignment effect [66]. If it jumps in the northeastern direction, then the chance of killing its prey is 82%. 
However, if it jumps in the opposite direction, then the chance of catching the prey is 18%. Consequently, it can be 
said that there are two to catch the prey. Fig. 1 shows the hunting behavior of the red fox. 
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X(it+1)=Dist_Fox_Prey(it)*Jump(it)*0.18 Wave sound of prey 


Dist Fox Prey(it)-Dist S T(it)*0.5 Jump(it)=0.5*9.81*t42 
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X(it) 
BN it+1)=Best_X+trand(1, dim)*MinT*a 
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\ Random walk J l J 


Exploration Exploitation 


Fig. 1 Red Fox Hunting Behavior: Exploration and Exploitation 


2.3 The FOX optimization algorithm 

In the beginning, FOX initializes the population, which is called the X matrix. An X is the position of red foxes. Then, 
the fitness of each search agent is calculated by using standard benchmark functions in each iteration. To investigate 
the BestFitness and best position (BestX), the fitness value of each search agent (each row in an X matrix) is compared 
to the fitness of other agents (other rows). BestFitness and BestX are conducted using a condition that is used to 
compare the fitness of the new row (f'itness;,4) and fitness of the previous row(fitness;) throughout iterations are 
returned. 


Then, to balance the exploration and exploitation phases, a condition is used with a random variable. This variable 
aims to divide the phases equally in terms of the number of iterations. Using a random variable called r, we assign 
50% probably to do either exploration or exploitation in FOX. In other words, nearly half of the iterations are used 
for exploration and the other half for exploitation. The proposed algorithm is essential in terms of balancing and 
avoiding local optima. Thus, a condition statement is used to divide the iterations equally for both exploration and 
exploitation. The variable a is used to decrease the search performance according to the BestX, and after each iteration, 
the value is decreased, which means that the agent pursues the prey better in every iteration. In addition to providing 
the condition to update the position, the fitness value affects the search agents to avoid local optima, because if the 
new position does not change, the exploration phase deactivates to allow other phases to be activated. The following 
two sections explain exploitation and exploration. 


2.3.] Exploitation 

In the exploitation phase, we have a condition relating to the chance of killing the prey as mentioned in Section 2.2. 
The value of the random variable p is in the range of [0,1]. Therefore, if the random number p is greater than 0.18, the 
new position of the red fox needs to be found. To find a new position, the distance sound travels Dist S T. the 
distance of the red fox from the prey Dist Fox Prev and jumping value Jumpj, must be calculated. Consequently, 
we generate a random number ranging from 0 to 1 for the sound travel time Time S T;,. The distance of the sound 
from the red fox is found by multiplying the speed of sound in the air Sp S with time sound travels time Time, S Tj, 
[67]: 


Dist S T, = Sp S * Time S Ti, (1) 


where the speed of sound in the medium Sp. S is equal to 343 in the air and Time. S Tj, is a random number in the 
range [0, 1]. Then, it is the number of iterations, which ranges from | to 500. However, another equation is constructed 
to find Sp S, which is based on the best position found so far by dividing the time sound takes to travel between the 
fox and prey. Time, S Tj, is a random and BestPosition;, is the best search agent in the population. Equation (2) 
shows how to find the speed of the sound Sp. S based on the best position. 
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BestPositionj, 
Sp S = ——————— 


TU Time Sh, @) 


To find the distance sound travels, equation (1) is used, which means that distance (Dist S Tj,)- velocity (Sp. S) * 
time (Time, S T;,). Consequently, the distance of the fox from the prey; Dist Fox Prey;, can be calculated by halving 
Dist S Tj. Therefore, to calculate the distance between sensor and object in physics, the distance of sound travels is 
divided by 2 since the distance between the sensor and the object is one-half of the distance that has been found by 
the sound wave [67]. Fig. 2 presents the sensor, which sends a sound wave signal to an object and then receives the 
signal. Therefore, the sensor multiplies the distance of the sound travels time by 0.5 or 1⁄2. Both multiplying by 0.5 or 
dividing by 2 are the same because half of the sound wave is considered [67]. So equation (3), is taken from techniques 
of sensors and objects to find the distance from Fig. 2: 


Dist Fox Prev = Dist S Tj, * 0.5 (3) 


Original wave 


Sender 


—-ODz3 7. Receiver 


Reflected wave 


Fig. 2 Calculating distance by ultrasound sensor [67] 


After finding the distance between a fox and the prey, the red fox needs to find a new position so that the red fox 
requires to jump to catch the prey. Therefore, the fox needs to calculate the jump height Jump,,. Thus, Jump;, can be 
calculated by the following equation: 


(4) 
Jump, = 0.5 * 9.81 « t? 


where 9.81 is the acceleration due to gravity and f is equal to the average time that sound takes to travel and it is 
squared because of the up and down steps in the jump. The time transition ¢t value is calculated by dividing the sum 
of the Time S Ti, to dimensions. Equation (7) shows the tt and MinT calculation. Then average time t is found by 
dividing tt by 2. The average time and gravity are multiplied by 0.5 because the jump value takes two different times 
to go up and down. As a result, both gravity and average time must be multiplied by 0.5. Then, the Jump value is 
multiplied by Dist Fox Prev and c,. The variable c,has a range, which is within [0, 0.18] when the red fox jumps 
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to the northeast direction. The equation below shows the calculation of the new position of the red fox if the p, which 
random number in the range [0,1], is greater than 0.18. 


X(it+1)= Dist_Fox_Preyj,* Jumpi, * cy (5) 


Both Equations (5) and (6) are used to find a new location for the red fox. Only one of them is executed in each 
iteration because of using the p condition. However, the only difference is in the second part of the p condition in 
Equation (5). If the else condition of the p is different, the equation (5) is multiplied by c; instead of c, if the value of 
p is less than or equal to 0.18. Therefore, if the p is greater than 0.18, then the new position is calculated by equation 
(5). However, if the value is less than 0.18, then the new position is calculated by equation (6). The range of c; is 
within [0.19, 1]. 


Xit+1) = Dist Fox Prey * Jumpit * c; (6) 


The value of c, and c; are 0.18 and 0.82 respectively. These values are based on the jump movement of a red fox, 
which either jumps to the northeast or the opposite. Therefore, if the p value is greater than 0.18, it means that the red 
fox jumps to the northeast direction. Consequently, to find a new position both Dist Fox Preyj, and Jumpi, are 
multiplied by c4. Accordingly, chances of exploiting a new position are great and the red fox goes toward global 
optima. However, the red fox jumps in the opposite northeast direction, if the p value is less than 0.18, this means that 
the chances of killing prey are low (%18). Hence, both Dist Fox Prev and Jumpy are multiplied by c;. 


2.3.2 Exploration 

To control the random walk, the fox searches randomly in this phase according to the best position of the fox that has 
been found so far. The fox does not have a jumping technique in this phase because it has to walk randomly to explore 
prey in the search area. To ensure that the fox walks randomly toward the best position, a minimum time variable 
MinT and the variable a are used to control the search. Equations (7) and (8) show the calculation of the MinT and 
variables. MinT is calculated by finding the minimum of tt. 


sum (Times, (i,: )) 
tt = —————, MinT = Min(tt 
dimension em PAD (7) 


Summation of Time_S_T;, (i,:) is divided by the dimension of the problem to find the minimum time average tt. 


a RT Go) (8) 


where Max; , is the maximum iterations. Calculating both MinT and a variable has a vital effect on the search phase 
to move toward the solution that is close to the best solution. Using rand(1, dimension) to ensure that the fox walks 
stochastically to explore the prey. However, to strengthen the searchability of FOX, both MinT and a variable are 
used. Variable r, which is a random number, is also used to balance the exploration and exploitation phases. The best 
solution BestX;, that has been found so far has a great impact on the exploration phase. Equation (9) shows the 
exploration technique of the fox in searching for a new position in the search space X (j¢4.1). The equations in this phase 
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can be adapted to existing algorithms to enhance their performance. They can also be used to propose new 
metaheuristic algorithms. 


Xit+1) = BestXj, * rand (1, dimension) * MinT «a (9) 


The equations in both phases do not need any modification rather than adaptation to a specific problem while FOX is 
used to solve the multidimensional space problem. The details of the FOX can be seen in Algorithm 1 and Figure 3. 


Algorithm 1 FOX optimization algorithm 

1: Initialize the red fox population X i (i—1, 2, ....., n) 

2: While it«Maxit 

3: Initialize Dist S T, Sp S, Time S T, BestX, Dist Fox Prey, Jump, MinT, a, BestFitness. 


4: Calculate the fitness of each search agent 

5: Select BestX and BestFitness among the fox population (X) in each iteration. 
6: Ifl fitness;^ fitness;,, 

Te BestFitness=fitness;.4 

8: BestX=X(i, :) 

9: Endif1 

10: If2 r >= 0.5 

11: If3 p >0.18 

12: Initialize time randomly; 

13: Calculate Distance_Sound_travels using Eq. (1) 
14: Calculate Sp_S from Eq. (2) 

15: Calculate distance from fox to prey using Eq. (3) 
16: Tt= average time; 

IE T=Tt/2; 

18: Calculate jump using Eq. (4) 

19: Find X (j¢41) using Eq. (5); 

20: Elseif p <=0.18 

2s Initialize time randomly; 

22: Calculate Distance Sound travels using Eq. (1) 
235 Calculate Sp_S from Eq. (2) 

24: Calculate distance from fox to prey using Eq. (3) 
E Tt= average time; 

26: T=Tt/2; 

DLT Calculate jump using Eq. (4) 

28: Find X (p41) using Eq. (6); 

29: EndIf3 

30: else 

31: Find MinT using Eq. (7) 

32: Explore X (41) using Eq. (9) 

33: 

34: EndIf2 

35: Check and amend the position if it goes beyond the limits 

36: Evaluate search agents by their fitness 

3/8 Update BestX 

38: it=it+] 


39: End while 
AO: return BestX & BestFitness 
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Initialize Population 
Calculate fitness of search agent 


Search for new X using Eq. 9 Find new X using Eq. 6 


Amend X if they are beyond 
the limitation 
Calculate fitness and 
return the best fitness 


Find new X using Eq. 5 


Fig. 3 Flowchart of FOX 


Overall, FOX begins by initializing the population of the red fox randomly. Then, in the first iteration, the population 
is checked whether the position of each red fox is inside the boundary of the benchmark function or not. After that, 
the fitness value of the benchmark function is calculated based on the row of the population. Following these steps, 
the BestFitness value and the best position (BestX) are selected. Then, a condition is started by comparing the random 
number r. If it is greater than or equal to 0.5, then the exploitation phase is activated. Also, a p condition exists inside 
the exploitation phase. If the p value is greater than 0.18, then the new position of the red fox is found based on 
Equations (1), (2), (3), (4), and (5). However, if p value is less than or equal to 0.18, then the new position is calculated 
based on Equations (1), (2), (3), (4), and (6). If r is less than 0.5 in the else condition, the exploration phase is activated. 
Therefore, the new position is found based on the best position, random number, and multiplying MinT variable by 
variable a. Thus, the BestFitness is returned in the first iteration. After modifying the population, the same steps are 
done repeatedly in the second iteration to find the best fitness and best position. 


Regarding the computational complexity of FOX: each iteration, has a time complexity of O(SearchAgents * D * it) 
where SearchAgents is the population size, D is the dimension of the problem, and it is the number of iterations. 
Therefore, it can be said that FOX has an O(n?) as time complexity. Also, the space complexity of FOX is calculated 
based on the vectors and matrices that can be seen in Algorithm 1. So, FOX has an O(n’) space complexity for each 
iteration. 
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3. Results and Discussion 

To ensure that the proposed algorithm works well, it must be tested by using different benchmark test functions. Five 
classical benchmark functions and the recent benchmark test functions are used to evaluate FOX performance. 
Besides, the results of the proposed algorithm are compared to Butterfly Optimization Algorithm (BOA), DA, ChOA, 
FDO, WOA, and GA. Then, these results are statistically compared to each other to determine whether the result is 
significant. Therefore, to ensure that the FOX works well to solve real-world applications, the FOX is used to solve 
the pressure vessel design problem and economic load dispatch problem. 


3.1 Classical benchmark functions 

Five different classical benchmarks were chosen to evaluate the performance of FOX [68]. These benchmark 
functions are divided into unimodal and multimodal functions. Unimodal functions are better for evaluating the 
exploitation phase because they have one local optimum. However, multimodal functions are used for the exploration 
phase due to having more local optima [69]. Table 1 illustrates both unimodal and multimodal functions. 


Table 1 Five Classical Benchmark Functions. 


No. Function Dimension Search fmin 
range 
Unimodal Benchmark Functions 
n 
1 fi) = 2: x; 30 [-100, 100] 0 
i-1 
n-1 
2 [ROEY Huna xD" G- 1Y] 30 [-30, 30] 0 
i-1 


Multimodal Benchmark functions 


S l 
3 h=), -xsin( lal) 30 | £900, 500] | 418 9829x5 


n 
fo(x) = D [x? — 10 cos(2mx;) + 10] [-5.12 
4 i=1 30 ae 0 
5.12] 
5 fuo) : z ? pi (=) +1 30 [-600, 600] 0 
x) = —— x= cos | — -600, 
4000 Zu; | i1 Wi 


FOX algorithms were executed 40 times with 100 iterations and then compared to five different algorithms, such as 
DA, FDO, and WOA. Table 2 presents results that are conducted from papers [13], [14] and the achieving results of 
FOX. As a result, FOX obtained better exploitation ability in the first unimodal function while WOA achieved 
optimum results in the second unimodal function. FOX also outperformed well in functions 3 and 4. However, FDO 
achieved the third rank which has a superior result in only function 5. Regarding the three multimodal functions, FOX 
also outperformed the other algorithms in two functions and it means that FOX has greater exploration ability because 
of dividing nearly half of the iteration for the exploration phase and the other half for the exploitation phase. 
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Table 2 Comparison Results of Five Classical Benchmark Functions for FOX. 


BOA DA ChOA FDO WOA FOX 
F 
Avera Std Aver Std Avera Std Aver Std Avera Std Aver Std 
ge age ge age ge age 
Fi 1.01E | 1.66E | 2.85E | 7.16E | 6.8573 | 0.0000 | 7.47E | 7.26E | 1.41E | 4.91E 0 0 
-11 -12 -18 -18 E-49 0003 -21 -19 — 30 — 30 
7.600 | 6.786 | 27.154 | 0.0016 | 23.50 | 59.78 | 0.0725 | 0.3974 | 38.43 | 0.082 
PZ 180909 | 0.0215 5 4 6 241 10 83 81 7 37 471 
16.01 | 9.479 | 5.6843 | 0.0012 | 14.56 | 5.202 
F3 | 28.68 | 20.178 88 1 E-14 031 54 > 0 0 0 0 
F4 1.35E | 6.27E | 0.193 | 0.073 0 0 0.568 | 0.104 | 0.0002 | 0.0015 0 0 
— 13 —14 3 4 7 2 89 86 
—3628. | 5.1249 - = 
F5 NA NA 2857 383.6 8022 2285 2066 5080 | 695.79 6097. 387.2 
4 207 84 .16 68 x 942 


3.2 CEC2019 benchmark test functions 


This modern benchmark suite includes 10 functions that were developed at the CEC conference[70]. These ten 
functions were invented for an annual optimization competition called “The 100-Digit Challenge”. The 
dimensionalities of these functions are different, and all of them are scalable. Functions CECO1 to CECO3 have 
different dimensionalities [14]. However, the other functions have the same dimensionality, which is [-100,100]. Table 
3 presents the CEC2019 functions. 


Benchmark Functions. 


Table 3 CEC2019 Benchmark Test Functions with Their Dimensions and Search Ranges: CEC2019 consists of Ten 


No. Functions Fř = Fix) D Search Range 
1 Storn’s Chebyshev Polynomial Fitting Problem 1 9 [-8192, 8192] 
2 Inverse Hilbert Matrix Problem 1 16 [-16384, 16384] 
3 Lennard-Joes Minimum Energy Cluster 1 18 [-4, 4] 

4 Rastrigin’s Function 1 10 [-100, 100] 

5 Griewangk’s Function 1 10 [-100, 100] 

6 Weierstrass Function 1 10 [-100, 100] 

7 Modified Schwefel ’s Function 1 10 [-100, 100] 

8 Expand Schaffer ’s F6 function 1 10 [-100, 100] 

9 Happy Cat Function 1 10 [-100, 100] 
10 Ackley Function 1 10 [-100, 100] 


The FOX is compared with GWO, ChoA, GA, and FDO [13]. These algorithms are selected because they are very 
competitive and are commonly used to solve benchmark functions and real-world problems. Therefore, the FOX is 
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implemented in MATLAB code, which can be found at https://github.com/Hardi-Mohammed/FOX, and 30 agents are 
used by each algorithm to search in the landscape. This search is performed using 500 iterations. The results are 
determined after running each algorithm 30 times, and then the average results are used for comparison. The results 
of other algorithms such as FDO, BOA, and ChoA are conducted from [13], [71], [19]. The algorithms that are used 
in this paper have used the following parameters as can be seen in Table 4. 


Table 4 Parameter Setting. 


Parameters Numbers 
Iterations 500 

Run of Algorithms | 30 

Search agents 30 

Dimension Based on Table 3 


Table 5 shows that the FOX outperforms all other algorithms on four out of 10 functions except FDO. However, the 
GA and GWO have the second rank among the six algorithms, both GA and GWO have the best results in only three 
functions. ChoA and BOA do not have great results against FOX, GA, and GWO. Therefore, ChOA and BOA have 
the worst result compared to other algorithms. 


In addition, FOX has superior results in 9 out of 10 functions compared to BOA, while it outperforms well in 8 
functions against ChOA. FOX also achieved the best results compared to GA in 7 functions. However, GWO has 
competitive results, so FOX is great in only five functions. 


FDO is a new metaheuristic algorithm and it achieved the best results. Therefore, FOX was also compared with FDO 
in Table 5, results proved that FOX obtained well in only two functions while FDO achieved optimum results in 7 
functions. Thus, FOX can explore efficiently compared with other algorithms and avoid local optima. It can also better 
exploit the best solution than the other algorithms in four functions. 


Table 5 Comparison of the GWO, ChOA, GA, BOA, FDO, and the FOX using 10 Functions from the CEC2019 
Functions. Results were conducted using 500 iterations with 30 runs. 


GWO ChOA GA BOA FDO FOX 
F 
Avg. STD Avg. STD Avg. STD Avg. STD Avg. STD Avg. STD 
1 2.13 3.07 4.24E+ | 9.67E+ | 5.32 7.04 5.89E+ | 1.14E+ | 1.66E+ | 1.54E+ | 2.58E+ | 22491. 
E+08 E+08 09 09 E+04 E+04 04 04 09 09 04 27 
2 18.343 | 0.0003 | 18.408 | 0.0185 | 17.35 | 17.349 | 1.89E+ | 2.91E- | 17.342 | 1.07E- | 18.344 | 0.0004 
9 04 3 8 02 01 01 9 14 2 25 
3 13.702 7.23 13.702 | 7.11E- | 12.70 | 13.702 | 1.37E+ | 6.17E- | 12.702 | 5.33E- | 13.702 7.11 
4 E-15 4 06 24 01 04 4 15 5 E-15 
4 300.98 | 686.81 | 5932.6 | 28552 6.23 61986. | 2.09E+ | 7.71E+ | 44.489 1.06E+ | 835.24 
15 53 2 07 E404 04 03 9 21.525 03 78 
0.2516 7.539 6.18E+ | 7.08E- 0.1669 0.8000 
5 | 2.4313 07 4.2094 | 0.8873 6 7.2765 00 01 1.1916 77 5.3148 19 
11.939 | 0.7307 | 12.154 7.400 1.18E+ | 7.71E- | 10.703 | 1.0345 1.4994 
6 4 45 4 0.6826 5 6.6877 01 01 3 59 5.0325 43 
7 534.76 | 292.02 | 1007.1 | 179.01 | 791.7 | 697.89 | 1.04E+ | 2.15E+ | 294.04 | 171.95 | 306.79 | 141.85 
55 04 3 66 42 03 02 8 19 4 88 
0.9939 6.100 6.34E+ | 3.59E- 0.6685 0.3819 
8 | 5.4021 56 6.7846 | 0.1562 4 5.8228 00 01 5.1731 29 5.4621 07 
14.732 | 49.951 5.31 2.27E+ | 8.11E+ 0.1645 0.4408 
9 8 42 449.27 | 245.49 E403 E403 03 02 2.5578 43 3.7959 76 
1 | 21.497 | 0.0685 | 21.498 | 0.0719 | 20.10 | 20.023 | 2.15E+ | 7.95E- 19.337 3.5913 | 20.981 | 0.0110 
0 3 13 5 5 59 01 02 À 28 8 96 
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3.3 Statistical analysis 

To evaluate the performance of the algorithm, statistical analysis is required to determine whether it is significant or 
not. Therefore, the Wilcoxon rank-sum test is used to check the p-value of the results. Table 6 shows that FOX is very 
competitive against all the algorithms. The FOX has a significant value in 7 functions against the GA and it has 9 
significant values against the BOA while it has 8 best results in terms of significance against the ChOA. However, it 
outperforms the GWO in only 6 functions. Overall, the FOX has superior results compared to the GWO, BOA, ChOA, 
and GA. These results prove that FOX is very competitive against common algorithms. 


Table 5 Statistical comparison (p-value) of the FOX against GWO, FDO, BOA, ChOA, and the GA using the Wilcoxon 
rank-sum test. 


F | FOX vs. GWO | FOX vs. FDO | FOX vs. BOA | FOX vs. ChOA | FOX vs. GA 
1 0.000 0.06 0.000 0.005 0.000 
2 0.008 0.102 0.000 0.003 N/A 
3 1.000 0.233 0.000 0.000 N/A 
4 N/A 0.008 0.000 0.004 0.000 
5 N/A 0.432 N/A N/A 0.000 
6 0.000 0.000 0.000 0.000 0.000 
7 0.002 0.06 0.000 0.000 0.000 
8 N/A 0.075 0.003 0.032 0.000 
9 0.000 0.62 0.000 0.001 0.000 
10 0.000 0.991 0.000 0.005 N/A 


3.4 Parameter sensitivity analysis and search mechanisms of FOX 

Both c4 and c; variables have a different range as described in the exploitation phase. c4 range value is [0, 0.18] 
while c; range value is [0.19, 1]. Therefore, to use a better threshold due to achieving better performance by FOX, six 
random numbers were chosen to evaluate the performance of FOX. The c4, random numbers are (0.16, 0.03, 0.16, 
0.02, 0.06, 0.18) and c; random numbers are (0.49, 0.65, 0.53, 0.37, 0.96, 0.82). By choosing each random number 
from both c, and c, respectively, the performance of FOX is evaluated. Fig. 4 presents that FOX obtained higher 
performance when the value of c, — 0.18 and c; — 0.82. FOX with both random values applied on all CEC2019 
functions. However, significant results can be seen in F4, F7, and F10. As a result, it can be said that both c4 and c, 
has a value of 0.18 and 0.82 respectively. Both c4 and cı shows as rand6 in Fig. 4 which shows that FOX achieved 
better performance with rand6. 
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Function Number 


Fig. 4 Evaluating Six Random Numbers 


Using the equations inside FOX and parameters guaranteed that the FOX achieves a better solution, as can be seen in 
Table 4, which includes five algorithms. Therefore, FOX is very competitive against the other algorithms. 


The variable and equations in FOX can be changed because the distance of the fox from the prey can be found by 
different equations. Additionally, the Jump variable can be found differently: Jump=0.5*9.81* t^2, which is used to 
find the jump value. However, it can also be used to subtract the current position from the best position. Therefore, by 
changing the jump equation, the results are changed. These parameters can be tuned as optimization parameters like 
a value, changing the equation for finding the distance of the fox from the prey. The value of a is one of the vital 
values that affect the exploration and exploitation phases. Therefore, by increasing the value, the result changes. 
Consequently, this variable can be used to tune the FOX. 


After implementing the FOX, three different versions of the FOX are implemented by changing the a value and jump 
equation. These modifications are based on the results in Table 4, which shows that FOX is not competitive against 
FDO. It means that FOX has limitations based on search capability and exploitation performance. Hence, the FOX 
exploration phase can be enhanced to be more competitive against other algorithms. The exploitation phase of FOX 
also can be improved based on changing the equations that are used to find the distance of red fox from prey and the 
jump equation. The convergence of the fitness value can be improved based on changing exploration and exploitation 
equations. As it can be seen from Fig. 6 FOX (a val+jump) version improves the FOX original. However, the results 
are required further improvements but it can be said that these parts can be considered as the limitation and 
improvement part for future work on the FOX. 


The FOX versions are the FOX with a val=4, the FOX using different jump equations, and both types of FOX 
combined (a val+jump equations). The results in the previous sections are related to the original FOX because the 
other versions do not improve against GWO, the WOA, CSO, and the GA. However, the original FOX can be 
enhanced based on the previously mentioned versions because they achieve better results than the original FOX. The 
results of these versions can be seen in Table 6. 


Table 6 Comparison of Four Different Versions of the FOX: FOX original, FOX (a val=4), FOX(jump), and FOX (a 
val+jump). 


Func. | FOX original | FOX a val=4 | FOX jump | FOX a val+jump 


1 2.58 E+04 3.17 E+04 2.12 E+04 1.53 E+04 
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2 18.3442 18.3442 18.3437 18.3443 
3 13.7025 13.7025 13.7024 13.7024 
4 1.06 E+03 189.7129 1.38 E+03 1.74 E+02 
5 6.295 3.4661 4.0394 3.4793 

6 5.0325 6.6608 5.4181 6.0872 

7 456.3214 371.0237 364.0703 423.5878 
8 5.6778 5.7067 5.7042 5.6309 

9 3.7959 3.3712 3.594 3.3702 
10 20.9878 20.9975 20.9994 20.9929 


To analyze the FOX versions in detail, two quantitative metrics were used. These two metrics are used to show the 
results as a graph to present the exploration and convergence rate of the search agents. Figures 5 and 6 show the results 
of these two metrics. 


The CECE2019 benchmark functions are also used in both metrics. Function 5 is used for the first metric, and then 
three functions were chosen for analysis randomly by the convergence metric, which are functions 1, 5, and 9. The 
FOX used 10 search agents using 100 iterations. Then, the results were obtained and put on a graph. 


The first experiment is used to illustrate the exploration capability of three different versions of the fox search agents 
and how they search through the search space. These versions are called a val, jump, and a val+jump. The version of 
the original FOX represented in this paper with a val has a value of this variable that is changed to 4. The second 
version is jump, which is the FOX with different jump functions. Finally, a val+jump is a combination of the previous 
two. In this technique, the best position of each iteration is recorded regarding its dimension. Therefore, the results of 
function 5 are as shown in Fig. 5. It can be seen that the foxes try to cover the landscape and then move toward 
optimality. Fig. 5 shows that FOX a val and FOX a val+jump achieve better results than FOX original and FOX jump. 
This shows the efficiency of search agents using a=4 because they search in a large area compared to the other types. 
Besides, the FOX and FOX jump do not explore efficiently compared to FOX a val. 


The second experiment is conducted by recording the fitness function of the three functions (functions 1, 5, and 9). 
Fig. 6 shows that the FOX versions start from a high fitness and then gradually decrease the fitness value until it 
reaches the optimum value. 
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Fig. 5 FOX Capability of exploration using Functions 5 
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Fig. 6 also shows that all four different versions of the FOX dramatically decrease and then stabilize during the 
iterations. However, the a val+jump version of the FOX performs well and achieves minimum results against the other 
versions in F5. Furthermore, the a val+jump version of the FOX has higher optimum results in F9, and the second- 
best version is the original FOX. 


Comparison of fitness F1 on differnt versions of FOX 
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Comparison of fitness F5 on differnt versions of FOX 
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Comparison of fitness F9 on differnt versions of FOX 
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Fig. 6 Convergence of FOX versions 


Overall, these improvements in the exploration capability and fitness value are related to the accuracy of the detecting 
techniques in finding the distance between the fox and the prey. Therefore, a better position is detected for the next 
iteration. As a result, the FOX is capable of searching in a landscape and can avoid local optima with efficient 
convergence toward optimality. FOX is better than the compared algorithms because it obtained a better average value 
when tested on five classical benchmark functions and CEC2019 benchmark functions. FOX is statistically effective 
because it has more significant values compared to the other algorithms. This algorithm is interesting because it has 
two different ways of exploitation based on the jump direction. It is also interested because FOX is inspired by a new 
idea which is the hunting behavior of the red fox. 


4. Engineering Case Studies to test FOX: 
In this section, two real-world applications are presented to be optimized by using FOX. The problems are pressure 
vessel design and economic load dispatch problems. 


4.] Solving a Pressure Vessel Design Problem 

Pressure vessel design is a classical engineering problem that has been solved by researchers using different algorithms 
[1], [72], [73]. Finding the optimal solution is the aim of this problem for three aspects of a cylindrical pressure vessel. 
These aspects are material, welding, and forming, and they should be optimized. This problem has four variables that 
can be optimized: the shell thickness T,, inner radius R, cylindrical length section without counting the head L and 
head thickness T,. Figure 7 illustrates the variables of the problem that can be optimized. However, the four 
constraints must be tested before optimizing the variables. The following shows the problem equation and constraints. 
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Fig. 7 Pressure vessel design problem [74] 


n = 1,2,3,4 


X = [x1x2x3x4] = [TsTh R L], 


f (3) = 0.6224x xax, + 1.7781x x2 + 3.1661x2x, + 19.84x2x,, (10) 


The variable limits are 0 < x <99, 


10 € x, < 200, 


10 € x, < 200, 


These are subject to the following constraints: 


gi (2) = -x + 0.0193x, < 0 (11) 

92(X%) = —x4 + 0.00954x, <0 (12) 
93(X) = —mx3x, — sexs + 1,296,000 <0 (13) 
ga(Z) = x4 +240 <0 (14) 


The WOA, WOAGWO, GWO, and FDO were used to solve this problem [51]. To prove that FOX has better 
performance against other algorithms, FOX is used to solve this problem. Table 7 shows that the FOX has high 
efficiency against the other algorithms. Overall, solving this problem is another proof that FOX is very competitive. 
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Table 7 Solving the Pressure Vessel Design Problem 


Run No. WOA WOAGWO FDO FOX 
1 13382.3 16697.97 377819.2 6866.846 
2 8106.258 13691.68 101800.8 7062.53 
3 13210.93 11391.75 57267.63 6535.226 
4 8369.843 11718.12 7556.644 7390.026 
5 8398.161 9536.122 6873.656 9511.966 
6 7438.549 11639.14 6777.363 14667.82 
7 19753.58 11127.94 19710.33 20734.75 
8 8205.785 15935.19 38812.92 24518.77 
9 8439.767 13165.6 8933.924 6346.123 
10 7966.724 10444.29 18258.82 6604.371 
11 8404.567 11303.79 116160.9 8786.426 
12 11396.46 16382.91 51808.24 17404.11 
13 10763.95 10050.45 83474.28 10341.99 
14 7714.56 11063.28 7538.922 8306.025 
15 159506.7 10520.38 47368.41 25313.13 
Avg. 20070.54 12311.24 63344.13 12026.01 
Std. 37401.74 2261.521 90960.63 6544.212 


4.2 Solving Economic Load Dispatch Problem 

The economic load dispatch problem is an optimization problem in electrical power. The main aim of optimizing this 
problem is to minimize the energy production cost while considering the load demand within various equality and 
inequality constraints [75]. Therefore, minimum fuel should be used to optimize the power generation unit to reduce 
the operation cost of generating energy [76]. Equation (15) shows the function of ELD. 


n 
C= 2,600 a 


Where C, is the total cost of fuel, P; is the actual power generation by generator i, C; is the required cost by generator 
and n is generator numbers. Therefore, the quadratic function used to represent C; as follow and the equation (16) has 
to be optimized: 


n 
i= 2. a;P? + b;P; + Ci (16) 
i=1 


= 


where each of aj, bj, and c; are used as coefficient cost for generator i. The above equation works under two 
constraints: the sum of the total power generator must be equal to the power demand with the power loss and the 
power generator must be inside the limitation of the power generator [77], [78]. 
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To solve the economic load dispatch problem with three power generators with a power demand of 150 MW. FOX 
was applied to solve this problem by using 1000 iterations with 50 independent runs. Table 9 shows the data by three 
different generators. Then, FOX results were compared against GWO, PSO, WOA, and FDO. All these algorithms 
are implemented to compare them against the FOX algorithm. Results proved that FOX achieved an optimum result 
in terms of generating powers while GWO, WOA, FDO, and FOX used nearly the same cost. Table 10 presents that 
FOX can obtain the required power at the minimum cost compared to other algorithms. 


Table 8 Power Limits and Cost Coefficients for Three Generators ELD. 


Generating Unit | Pmin | Pmax | a($/MW2) | b($/MW2) | c($/MW2) 
P1 10 85 0.008 7 200 
P2 10 80 0.009 6.3 180 
P3 10 70 0.007 6.8 140 


Table 9 Comparison of FOX with other Algorithms for Three Generators ELD Problem when PD = 150 MW. 


Unit/Algorithm GWO PSO WOA FDO FOX 

P1 (MW) 31.94 60.0345 | 31.938 32.665 | 31.937 
P2 (MW) 67.284 25.6626 | 67.284 65.489 | 67.277 
P3 (MW) 50.777 67.2313 | 50.778 51.846 | 50.785 
Total Generated Power (MW) | 150.001 | 152.9285 | 150 149.999 | 152.6089 
Cost ($/hr) 1579.698 | 1637.084 | 1579.699 | 1579.87 | 1579.699 


5. Result analysis 

The FOX algorithm demonstrated better exploration and exploitation when solving five classical benchmark functions 
compared to BOA, DA, ChOA, and WOA. It can be said that FOX has a great characteristic in terms of exploitation 
because it has optimum results compared to other algorithms in unimodal functions. Table 2 presents the exploitation 
capability of FOX for solving unimodal functions. FOX also benefits from superior exploratory behavior because it 
achieved the best results in solving multimodal functions in Tables 2 and 5. Therefore, the reason behind this efficiency 
is that the FOX avoids local optima because it uses variable a that affects the exploration phase to go toward the best 
solution. Dividing the iterations so that nearly a half is used for exploration and the other half is used for exploitation 
by using a condition with a random variable has a great impact on achieving a better balance between the two phases. 
Having an exploration phase based on best position, variable a, and MinT variables would have a great impact to 
explore solutions that are going forward with the global optima. Furthermore, finding the distance between fox and 
prey with a fox jump method is a vital method that is used to catch the prey in the exploitation phase. As a result, FOX 
achieved better performance against FDO, BOA, GA, GWO, and ChOA. 


FOX optimizer uses a physical method to find its distance from prey by retrieving the sound of the prey under the 
snow and then trying to use a jump mechanism to catch prey. Finding distance from prey and jump method has benefits 
over LOA, Tigers, and GWO because these algorithms use social hierarchy and group hunting mechanisms. Despite 
these technics, the FOX optimizer achieved better performance compared to GWO[11], Tigers [39], LOA [38]in 
solving benchmark functions. 


Three different versions of FOX were compared in Table 6 to show that FOX can be improved by changing the value 
of variable a and jump equation. These versions can improve FOX but they are not statistically significant. Therefore, 
it can be used in terms of modifications of FOX. The reason is that the equations of jump can be identified by different 
methods in the future and also the value of a variable has a great impact on the exploration of FOX, so by increasing 
the value of a, the performance of FOX is improved. 
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FOX is also used to solve pressure vessel design problems and economic ELD problems. Then, it achieved better 
performance against WOA, WOAGWO, GWO, PSO, and FDO. However, FDO performance is better than FOX in 
solving CEC2019 benchmark functions. These results prove that an algorithm can work well to solve problem A while 
it does not work well for problem B. The reason is that FOX has better exploration capability to avoid local optima 
and divide the iterations into both exploration and exploitation phases. 


Overall, it can be said that FOX can avoid local optima and it also has proper balancing between exploration and 
exploitation by using random variable r. FOX also mimics the behavior of the red fox by using physical technics to 
create different equations from the red fox behavior. Finding the distance between the red fox and prey by taking 
advantage of the physical experiment in Fig. 2 improves the performance of FOX. 


6. Conclusion 

To conclude, this paper presented a new metaheuristic algorithm inspired by foxes’ hunting behavior in nature called 
FOX. It uses a random walk based on the best position in the exploration phase, while it uses a distance measurement 
method to find the distance between fox and prey. It also uses jump action to catch the prey in the exploitation phase. 


The FOX was tested by using five classical benchmark functions and CEC2019 benchmark functions, then compared 
to the GWO, ChOA, GA, PSO, BOA, and GA. Then, statistical results showed that FOX has the first rank and is better 
than the compared algorithms. The reason behind the performance of FOX is that FOX uses jump technics to catch 
prey. The other reason is that FOX searches based on best position, a variable, and MinT. 


To provide parameter sensitivity analysis and discover more insights into the search mechanisms of FOX, exploration 
metric and convergence capability metric were used to present the performance of FOX. Three different versions of 
FOX are used by the metrics and the results showed that the FOX a val and FOX (a val+jump) have a great ability to 
explore new positions against the FOX original. As a result, both of these versions perform well in terms of 
exploration. The convergence metric is implemented by using Functions 1, 5, and 9. Jump-a val FOX versions showed 
excellent convergence and achieved better results than the other versions in functions 5 and 9. Therefore, FOX requires 
further research improvement. 


The FOX was also applied to solve a pressure vessel design problem and economic load dispatch problem. The FOX 
achieved results better than WOA, WOAGWO, GWO, KMGWO, and FDO. FOX generated optimum power results 
because it has efficient exploration and exploitation capabilities. The reason behind the high performance of the FOX 
is that the FOX uses a proper balance technique between the exploration and exploitation phases to avoid being trapped 
in local optima. Decreasing the value of a specific variable throughout the iterations enhanced the exploitation phase. 


This algorithm can be further extended since the FOX works for single-objective optimization; multi-objective 
optimization will be studied in the future. The FOX can also be used to solve optimization problems in planning and 
medical problems. 
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